## fully independent model

model {
	for (i in 1:(H+1)) {
		y[i] ~ dbinom(psi[i], n[i])
	
		psi[i] ~ dbeta(0.5, 0.5)
	}
}